Consequences of Anomalous Diffusion in Disordered Systems Under Cyclic Forcing 
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We use numerical simulations to study the behavior of 2D frictionless disk systems under cyclic 
shear as a function of reversal amplitude y r . Our studies focus on mean bulk and disk dynamics. 
These measurements suggest a crossover from a subdiffusive, 7 r dependent regime to a regime where 
the grain motions are diffusive, with properties dependent only on total shear strain. We discuss 
model stochastic processes that are consistent with these observations. Finally, we introduce a modi- 
fied Mean-Squared Displacement (mMSD) which takes into account the motion of the neighborhood 
of nearby grains and yields new insights into local displacement fluctuations. We find that scaling 
properties of the displacement distributions are consistent with well studied stochastic models of 
anomalous diffusion and suggest scale-invariant cage dynamics. 
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Introduction — In recent years there has been a good 
deal of interest in amorphous systems under oscillatory 
driving [THE]. These works focus on the property of re- 
versibility of particle dynamics, or the dependence of dy- 
namical behaviors on the distance to tf>j, which is the 
density at which the packing becomes "jammed." How- 
ever, what has been less well studied is the systematic 
dependence of bulk and particle dynamics of amorphous 
systems on driving amplitude, for systems that are not 
strictly reversible. In this study, we simulate an amor- 
phous system under slow, cyclic shear for a range of rever- 
sal amplitudes 7 r (the strain at which shear is reversed). 
We find that for our simple model, a crossover exists in 
the behavior of bulk and particle dynamics from a regime 
that is dependent on j r to one that is not. This crossover 
appears to influence a wide range of properties, including 
shear stress, packing fraction, and particle displacements. 

Finally, we address the microscopic origins of the j r de- 
pendent diffusive properties. We find a j r dependence in 
the exponent /3 of the mean-squared displacement (MSD) 
(r 2 ) (x t@ . To limit the possible underlying stochas- 
tic processes that can potentially explain the observed 
sub-diffusivity, we analyze the the scaling properties of 
probability distribution functions (PDFs) as a function 
of distance and time for all grains. We observe that the 
exponent j3 exclusively dictates the 7 r dependence of the 
PDFs, as well as the MSDs and the dependence of bulk 
properties on j r . 

Methods — We study a minimal model system with the 
use of molecular dynamics (MD) simulation. Bidisperse 
(1:1) frictionless disks, with larger disks having 1.4 times 
larger diameter than the smaller, interact via a linear 
spring potential and damping[7] that is dependent on the 
difference in velocity between contacting grains. Lees- 
Edwards boundary conditions are used to apply a shear 
strain of 7. It is important to note that we choose not 
to include a coupling to an affine displacement field since 
we only intend to capture interactions between grains 
at contact, and not interactions with a constituent fluid. 
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FIG. 1: The schematic diagram illustrates the shear as applied 
to the 2D model granular system. The initial state has a 
repeat cell which is square with linear size L. The strain is 
applied by shifting the the lower and upper adjacent repeat 
cells horizontally by up to Ax r as per Lees-Edwards boundary 
conditions. Then, the shear strain Ax r is reversed back to 
zero. 

Packings are kept at approximately fixed pressure during 
shear by uniformly modulating the grain sizes in response 
to changes in pressure. At each simulation time step, 
the change in pressure is calculated, and if the pressure 
decreases, then grain diameters are increased, otherwise 
the grain diameters are decreased (the pressure is kept 
approximately fixed to the same value throughout this 
work). Under these conditions, we study 8000 grain sys- 
tems over a range of j r = ^jf- from 0.05 to 0.225. 

Results — We begin by presenting some of the dynam- 
ical properties of our system. In particular, we look at 
the mean squared displacements (MSD) of the grains. 
We will always report MSDs on a cycle- by- cycle basis, 
that is to say, the MSD is calculated at the end of each 
cycle, where for reversible dynamics the system should 
return to zero net displacement. We find that for all 7,. 
studied grains are at least slightly displaced relative to 
the beginning of the cycle, with increasing MSDs as a 
function of cycle number (where cycle number replaces 
time). First, we simply report the MSD relative to the 
initial (isotropic) state. These MSDs are noisy, making it 
difficult to extract a single j3 characterizing the evolution 
of the grain displacements. As a result, we also average 
over intervals of up to 200 cycles, each separated by 10 
cycles in the oscillatory shear trajectory, and produce a 
time (or cycle) averaged MSD (see figure [2]). These av- 
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FIG. 2: (Top) MSD measured relative to the initial, isotropic 
state. Different colors correspond to different reversal ampli- 
tudes, with blue being the lowest and red being the highest. 
Solid black lines with a slope of 1 are added for reference. 
(Bottom, inset) Time (or cycle) averaged MSDs. The best 
linear fits are shown for each reversal amplitude. (Bottom, 
main) the slopes from the linear fits are shown as a function 
of reversal amplitude. 

eraged MSD curves exhibit a single exponent over the 
range of 200 cycles. Neglecting the first 10 cycles as the 
very short time dynamics appear to be somewhat distinct 
from the longer time scale dynamics, we apply linear fits 
to the MSD curves and extract slopes. These slopes as 
a function of 7 r provide the first evidence of a crossover 
between two regimes as 7 r is increased. For very small 
7 r , the dynamics are highly subdiffusive. However, for a 
7 r =0.15 the exponent f3 is nearly 1, suggesting diffusive 
dynamics. As j r continues to increase, j3 stays unchanged 
at approximately 1. 

This crossover in grain dynamics as 7 r is increased also 
manifests itself in the bulk properties of the shear flow 
(figure pjj). During each cycle, the packing fraction <f> in- 
creases to a peak value, which depends on y r , and then 
decreases again. Averaged over many cycles, this process 
of compaction followed by dilation appears to begin and 
end at the same value of <fi and for smaller "f r exhibits a j r 
dependent strain value at which the peak <j> is reached. 
This strain at peak is plotted against 7 r in figure |3j 
and a similar crossover from a 7 r dependent regime near 
7 r =0.15 is again observed. 
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FIG. 3: Packing fraction (top, inset) and shear stress a xy (bot- 
tom, inset) averaged over many cycles as a function of strain 
7, for varying 7,-. Colors range from blue (smallest y r ) to red 
(largest 7 r ). (Main) f r dependence for strain (top) at <f> and 
slope of a xy (bottom) at reversal. 

A second observed bulk property that shows this 
crossover behavior is the shear stress o~ xy . The rate of 
change of a xy again shows a 7 r dependence. We illustrate 
this result by calculating the difference in a xy between 
the last two average strain steps before reversal, which 
again shows a crossover from a 7 r dependent regime near 
7 r =0.15 (see figure |3|. We also note that averaged over 
many cycles, a xy is quite hysteretic. Very quickly the 
stress state of the system even at the beginning and end 
of each cycle is anisotropic, oscillating on average be- 
tween two extremal values of o~ xy . 

Finally, we note that long time transients exist partic- 
ularly for small j r . The essential conclusion from these 
results is that a crossover in dynamics and bulk proper- 
ties as a function of j r occurs. Our finding is insensitive 
to the transient: <p an d °xy increase for lower 7 r over 
many cycles initially, but the strain at peak <f> and slope 
of a xy quickly become independent of cycle number. As 
for the dynamics, the MSD curves shift to lower squared 
displacements, but the slopes do not change significantly. 

Connecting to Models of Subdiffusion — With such a 
strong correspondence between changes in bulk proper- 
ties, and a crossover from a subdiffusive regime to simple 
diffusion, our next goal is to identify models that can 
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capture the relevant subdiffusive dynamics. A concise 
survey of stochastic processes that result in subdiffusive 
dynamics can be found in [8]. These models are all able 
to exhibit a subdiffusive regime described by an MSD 
which is oc f 9 , but are motivated by different physics. 
In brief, there are two basic classes, fractional Brown- 
ian motion (fBm) and fractional time processes (ftp). 
Both exhibit self-similarity under time dilation. How- 
ever in the case of fBm a time dependent diffusion coef- 
ficient is introduced, while for ftp a heavy-tailed distri- 
bution of waiting times is introduced. As a result, the 
latter process is considered non-Markovian. In addition, 
both processes can be generalized to fully fractional forms 
by introducing heavy-tailed displacement distributions in 
analogy to Levy flights. The PDFs for all of these pro- 
cesses must obey a particular scaling form that reflects 
the self-affinity: 



t H P{r,t) oc ip (£ 



t H 



(1) 



The function tp (£ = r/i H ) will have different forms for 
the different models discussed. Importantly, the self- 
similarity exponent H — (3/2 in the cases of fBm and 
ftp, while the fully fractional processes do not generally 
satisfy these relationships. Motivated by this fact, we 
scale PDFs for each of our trajectories at different j r us- 
ing the exponents gained from the time averaged MSDs. 
As an example, Figure [4] shows PDFs for the time inter- 
vals At = 10, 50, 100 for the case of j r = 0.075, scaled 
according to equation [T] with H = (3/2. Collapses for 
other 7 r are of similar quality. We conclude then that the 
fully fractional processes are not the appropriate model 
for the subdiffusive dynamics observed here, but that ftp 
and fBm are both viable candidates. 

Modifying the MSD — To further differentiate between 
ftp and fBm, we consider the physical implications of 
waiting time distributions. Waiting times can be mod- 
eled by continuous time random walks (CTRWs) in which 
the random walker has the ability to wait in one place 
rather than moving. When the amount of time associated 
with waiting in one place has a heavy-tailed distribution, 
a ftp results. This description shares similarities with the 
'cage breaking' scheme used to describe long relaxation 
times in granular and glassy systems. In this scheme, 
particles such as the grains in our simulation 'rattle' in- 
side a cage of their neighbors on a shorter time scale, 
intermittently escaping out and into a new cage. The 
time spent in cages is akin to a 'waiting time.' The ftp 
model would suggest that these waiting times should be 
heavy-tail distributed if the cage effect is to result in sub- 
diffusive MSDs. However, the notion of a cage is difficult 
to define in a system where all of the particles are in 
motion, so that the cage itself is a transient feature. 

To further explore the validity of this cage breaking 
scheme, we modify the MSDs (see figure [5]). To high- 
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FIG. 4: (Top) The unsealed displacement distributions for 
cycle intervals At = 10, 50, 100 (for black, blue and red). 
(Bottom) the displacement distributions scale according to 
equation [I] using the slope from the time averaged MSD for 
j r = 0.075. 



light motion relative to the cage, we calculate a modified 
MSD (mMSD) for a given grain i with respect to a neigh- 
borhood of nearby grains. In the isotropic state of the 
packing (7 = for the first cycle) the 6 nearby grains 
are found for grain i. A centroid for these nearby grains 
is calculated at every strain step, and the displacement 
of grain i relative to the centroid is calculated. In this 
way, we characterize the motion of a grain relative to 
its nearby grains, which provides better insight into the 
nature of the local fluctuations in grain displacements. 
Figure [6] shows that the mMSD curves are very smooth, 
with total squared displacements almost an order of mag- 
nitude smaller than those found with the MSD. These 
results suggest that the motions of the grains are coordi- 
nated with the motions of their neighborhoods of grains, 
which is itself changing size and shape over time. Inter- 
estingly, unlike the MSD, the mMSD exhibits a contin- 
uous change from subdiffusive to superdiffusive motion, 
with superdiffusion for j r > 0.15. Understanding the ori- 
gin of superdiffusion in this context is a goal for future 
investigations. 

Discussion — We have presented a systematic study 
of reversal amplitude j r in a model amorphous system 
under cyclic shear. Results for bulk and grain dynamical 



FIG. 5: The above diagram illustrates the modification to the 
standard MSD approach. In the standard approach, grain 
displacements are calculated relative to their own locations 
at a previous time (the displacement vector, in red, is la- 
beled 'MSD'). For the modified MSD (mMSD), described by 
the inset, particle displacements are calculated relative to the 
centroid (green star) of the nearby 6 grains (in green). As 
the green grains move, the centroid moves as well, and the 
displacement of the blue grain is measured relative to this 
point. The displacement vector associated with the mMSD is 
labeled in the bottom left. 

properties suggest a crossover in behaviors from subdif- 
fusive to diffusive behavior near "/ r = 0.15. Interestingly 
the MSD exhibits a novel crossover in particle dynamics 
as a function of increasing j r , that is distinct from the 
crossover from 'caged dynamics' to diffusion as a func- 
tion of time reported in thermally driven systems such 
as colloids [9]. Temperature drives the system at the par- 
ticle scale, while cyclic shear drives the system at the 
boundary, i.e. the driving force is applied at large scales. 

Together these results indicate that driving from large 
scales rather than at the particle level fundamentally al- 
ters particle motion: the observed displacement distri- 
bution and its scaling behavior are not consistent with 
a simple cage-breaking model that is typically used for 
thermal systems. Furthermore, a mMSD measurement 
that characterizes grain motion relative to its neighbors, 
suggests that a 'cage breaking' scheme may not pro- 
vide the appropriate description of grain dynamics under 
cyclic shear: instead of rattling in their cage, grain MSDs 
are almost an order of magnitude smaller when measured 
relative to their cage. 

To expand on this idea further, we look at displace- 
ment distributions relative to the neighborhoods. The 
approach used here is to find the centroid of nearby grains 
during each cycle, and find the relative motion after one 
cycle. For an example trajectory with 7 r = 0.075, mod- 
ified displacement distributions and standard displace- 
ment distributions for At = 1 are shown in figure [7] All 
displacements are much smaller for the modified case, 
with a mean distance after one cycle approximately an 
order of magnitude smaller for the displacements relative 
to the neighborhoods. If we consider the amount of time 
a grain spent inside its neighborhood (defined as being 
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FIG. 6: The mMSD, relative to the packing configuration at 
t = 0, is shown. Solid black lines have a slope of 1. (Inset) 
MSDs for 7r = 0.225 for MSD (black) and mMSD (red). 

within the grain radius from the centroid) to be analo- 
gous to a cage lifetime, we find that the distribution of 
cage lifetimes appears to be scale invariant, with an ap- 
parent power-law distribution (inset of figure [7]) . This 
observation is consistent with the ftp model, where wait- 
ing times are assumed to be heavy tailed, although the 
exponent of the observed cage lifetimes (which appears 
to be 1) is not consistent with ftp. Still, the observa- 
tion that the cage lifetimes lacks a single characteristic 
time scale is at odds with the standard caging scenario, 
such as presented in [5], where the cages have a single 
characteristic lifetime for a particular density. 

Similar questions concerning the origins of subdiffusion 
in dense systems are ubiquitous in biological physics as 
well. Subdiffusion has been observed in for instance, the 
cytoplasm of living cells [TOT ITT] and mRNA dynamics in 
E. Coli [T5]. Interestingly, a discussion has emerged as 
to whether ftp or fBm is most appropriate in the case of 
mRNA dynamics as well [13] . 

Beyond the implications for granular dynamics sug- 
gested by the ftp model, such as long-time correlations, 
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FIG. 7: The single cycle displacement PDFs as a function 
of displaced distance r is shown for 7 r = 0.075, for both the 
displacements (black) and the relative displacements (red). 
(Inset) Waiting times, as defined in the text, appear to have 
a heavy tailed distribution (again the example of 7 r = 0.075 
is shown). Dashed line has a slope of -1. 
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the observation that the crossover in the slopes of the 
MSDs from less than 1 to 1 is consistent with the 
crossover in bulk density and stress properties has futher 
implications. Since the origin of the collapse of the dis- 
placement PDFs is the self-affinity of these PDFs, and 
the self-affinity exponent controls the slope of the MSDs, 
we suggest that this self-affinity may be the origin of the 
crossover in all bulk properties as a function of j r , al- 
though the reason for such a strong dependence on the 
self-affinity exponent H is not yet clear to us. We sug- 
gest that a more complex picture of particle dynamics 
in cyclic, athermal, flows is needed, which will require 
more focus on the importance of global boundary forcing 
versus particle scale thermal forcing as well as driving 
amplitude. 

The computational work reported here was done 
largely on the University of Maryland High Performance 
Computing Cluster (HPCC). This work was supported 
by NSF grant DMR0907146. 
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